Lattice-corrected strain-induced vector potentials in graphene 
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The electronic implications of strain in graphene can be captured at low energies by means of 
pseudo vector potentials which can give rise to pseudomagnetic fields. These strain-induced vector 
potentials arise from the local perturbation to the electronic hopping amplitudes in a tight-binding 
framework. Here we complete the standard description of the strain-induced vector potential, which 
accounts only for the hopping perturbation, with the explicit inclusion of the lattice deformations 
or, equivalent ly, the deformation of the Brillouin zone. These corrections are linear in strain and 
are different at each of the strained, inequivalent Dirac points, and hence are equally necessary to 
identify the precise magnitude of the vector potential. This effect can be relevant in scenarios of 
inhomogeneous strain profiles, where electronic motion depends on the amount of overlap among 
the local Fermi surfaces. In particular, it affects the pseudomagnetic field distribution induced by 
inhomogeneous strain configurations, and can lead to new opportunities in tailoring the optimal 
strain fields for certain desired functionalities. 

PACS numbers: 81.05.ue, 73.22.Pr, 71.70.Di, 85.35.-p 



I. INTRODUCTION 

Many surprising properties of graphene can, and have, 
been engineered to create novel physical scenarios in a 
condensed matter setting: the chiral and Dirac-like en- 
ergy dispersion near the Fermi energy leads to Klein 
tunneling 1 and Veselago lensing 2 ; relativistic Landau 
level quantization, observable even at room temperature, 
emerges under applied magnetic fields^. Graphene is also 
an extremely strong material, with an effective Young's 
modulus of 1 TPa, an intrinsic ultimate strength of 130 
GPa 4 , which has been seen to withstand elastic defor- 
mations up to 15- 20#l This unusually large range of 
elastic response, combined with graphene's intrinsic two- 
dimensionality and the peculiar coupling between elec- 
trons and deformations, make graphene an unparalleled 
condensed matter system in which to explore strain as a 
tool for tailoring electronic functionality. 

A dazzling glimpse of the feasibility and potential of 
strain-engineered graphene 6 7 has recently emerged with 
experiments reporting that certain strain profiles can in- 
duce Landau quantization and effective pseudomagnetic 
fields in excess of 300 T^E^. Such physics strongly en- 
courages the prospect of harnessing this unconventional 
interplay between graphene's unique electronic and im- 
pressive mechanical properties to control electronic trans- 
port in graphene devices^-4 

The key elements underlying these experimental ob- 
servations and the concept of strain-engineering can be 
readily captured from a tight-binding description of the 
electrons in graphene. Graphene has a honeycomb lat- 
tice (Fig. [IJa)) which can be decomposed into two tri- 
angular Bravais sublattices, the A and B sublattices, 



with nearest neighbor vectors of 5\ = a/2 x + \fZa/2 y, 
82 — a/2 x — \fZa/2 y, and £3 = —a x with a = 1.42 A. 
The two elementary lattice translation vectors (Ri) are 
also shown. The associated band structure was first 
calculated in 1947 by Wallace, considering only the tt 
band electrons^ 2 -'. It is characterized by two energy 
bands touching at the two non-equivalent corners (K 
and K')of the hexagonal Brillouin zone (BZ) 12 , shown 
in Fig. [lib). The positions of the Dirac points can be 
found by rotating K = V by 60 degree increments 

about the origin. The dispersion is conical near these 
points, mimicking the physics of 2D relativistic fermions. 
In undoped graphene, the Fermi energy lies exactly at the 
contact point between the two dispersing cones. At these 
so-called Dirac points, an effective low energy Hamilto- 
nian can be derived which coincides with the 2D Dirac 
Hamiltoniari^tt^] p Q r unstrained graphene, the three 
K points are equivalent, and likewise for the three K' 
points. 

After accounting for the effects of the distorted lat- 
tice, a tight-binding description is equally applicable 
to strained graphene. In Fig. [ljc), uniaxial strain in 
the armchair direction is shown. The distance between 
the nearest neighbors is modified, introducing a bond- 
dependent nearest neighbor hopping amplitude 16 . Ad- 
ditionally, the nearest neighbor vectors have been modi- 
fied as well. The essential consequences of these changes 
are the following^ (i) for any amount of strain, the 
Dirac points are displaced from the corners of the un- 
strained BZ and, furthermore, do not necessarily sit at 
the corners of the strained BZ; (ii) the gapless and con- 
ical nature of the energy dispersion remains robust, ex- 
cept when the deformation is so strong that the two in- 
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FIG. 1. (a) Unstrained graphene's real space lattice with la- 
beled nearest neighbor vectors (Si), labeled lattice translation 
vectors, Ri = na+ + ma_ with m and n integers, and black 
and gray dots representing the A and B sublattices, respec- 
tively, (b) The Brillouin zone of unstrained graphene with 
labeled high symmetry points, (c) The positions of the un- 
strained (black, solid) and strained (brown, dashed) nearest 
neighbors, and the corresponding Si and S[ for 20 % uniax- 
ial strain along the armchair direction, (d) The unstrained 
(black, solid) and strained (brown, dashed) Brillouin zone, 
also for 20 % armchair uniaxial strain, with labels for the 
now inequivalent Dirac points. 

equivalent Dirac points merge in a Lifshitz transition 
(but that probably requires strains of the order 20%, 
where the tight-binding description is not reliable any- 
more); (iii) at any finite density the Fermi line is de- 
formed from the isotropic circle to an elliptical shape, 
and two Fermi velocities can be defined along the princi- 
pal directions'^. All these modifications are significant 
and happen concurrently, and hence a complete descrip- 
tion of the electronic and transport properties of strained 
graphene requires their combined consideration. For ex- 
ample, the local shift of the Dirac point (i) can hinder 
or completely suppress electr onic propagation across re- 
gions of different strain state d 6 * 11 !. The anisotropy of the 
Fermi surface (iii) has direct bearing in measurable quan- 
tities such as t he ani sotropy in electrical resistivity^, op- 
tical absorptio n 18 * 21 ^, and the Raman signature of the 2D 
peajPBH. 

From the theoretical as well as technical point of view, 
the effects of strain are frequently considered indepen- 
dently, and one usually isolates the dominant effect for 
the physical observable of interest. Referring to the same 
examples above, the strain-induced corrections to optical 
absorption arising from inter-band transitions are insen- 
sitive to the absolute position of the Dirac point in the 



BZ, but strongly depend on the velocity anisotrop ^ 18 * 21 1 . 
Likewise, the dominant effect across a strain barrier will 
be the relative position of the Fermi surfaces in the two 
regions (since this essentially determines the phase space 
available for transmission), and in a first approach the 
anisotropy is usually neglected 6 -^, since the full descrip- 
tion would obfuscate the presentation of the problem. 

When the strain-induced shift of the Dirac points (i) 
is considered independently of (ii) a nd (iii) , it can be 
thought of as a pseudo vector potentia l 14 * 26 ! ^. This can 
be done because of the peculiar form of the strain correc- 
tions to the electronic dispersion in graphene. Electrons 
in strained graphene are still governed by a Dirac equa- 
tion, but one in which the strain modifications can be 
completely absorbed in the replacement p — ^ p — eA 
where A is the pseudo vector potential. This matches 
the conventional minimal coupling scenario, which means 
that the electrons respond to the deformation-induced 
perturbation as they would to an external magnetic field. 
The pseudo vector potential is related to the shift in the 
Dirac point, A/c^, through Akjj = — |A. This analogy 
between strain-induced and real magnetic fields means, 
for example, that the electronic energy levels can be 
quantized with a relativistic Landau spectrum just as if 
they were under a real magnetic field given by B = V x A 
(as long as this pseudomagnetic field is relatively constant 
on scales not smaller than the corresponding magnetic 
length)P. 

An omission in earlier work in the context of these 
pseudo vector potentials is the explicit consideration of 
the deformation of the lattice when computing the po- 
sition of the new Dirac points. Here we include this ef- 
fect, and show its importance in describing the absolute 
position of the Dirac points, and the resulting pseudo 
vector potentials. This yields new leading order terms 
in the strain-induced pseudo vector potential which are 
different at the three inequivalent Dirac points. Specifi- 
cally, we shall be interested below only in how strain af- 
fects the position of the Dirac point in reciprocal space. 
We also restrict the discussion to planar deformations, 
and hence i gnore effects that might arise in the presence 
of curvatur d 14 | 29 l We will detail the derivation of these 
terms in the Theory section, then, in the Discussion sec- 
tion, demonstrate their importance in describing the shift 
of the Dirac point in graphene, and, finally, illustrate how 
the inclusion of these terms affects the pseudomagnetic 
field distribution for particular strain profiles. 



II. THEORY 

Deformation of the lattice necessarily leads to modi- 
fied hoppings. If this is treated as a slowly varying per- 
turbation to the relaxed nearest-neighbor tight-binding 
Hamiltoniarf^, the Hamiltonian can be written as 

H = -t J2 a l b j ~ 3ti i a \ h i + Hc - w 

<i,j> <i,3> 
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Here the sum is over all nearest neighbor pairs, a\ is the 
creation operator for an electron on the A sublattice, bj 
is the annihilation operator for an electron on the B sub- 
lattice, t is the unstrained hopping amplitude, and Stij 
is the nearest neighbor, bond specific, change in hopping 
energy due to strain. The bond specificity of St can be 
made explicit by writing SUj = St(\Sij\). 

The importance of the lattice deformations is hidden 
in the phase factors that lead to the dispersion relation, 
and which arise upon writing eq. in Fourier space. 
To explicitly show this dependence, we work in the basis 
in which the phase factors in the Fourier transform are 
given by the positions of the strained sublattices^J. The 
creation/annihilation operators are written as 



(2a) 
(2b) 



The correct association of a finite pseudo vector poten- 
tial includes both the actual modification of the relative 
positions of the atoms, which enters in eq. (|2|, and the 
change in the hopping amplitudes, entering in eq. ([!]). 
Clearly, they are not controllable independently in an 
actual physical system, since the latter is a consequence 
of the former. In the past, the modification of the relative 
positions of the atoms has been ignored. This was done 
because in theoretical calculations one usually abstracts 
from the actual underlying lattice when considering the 
tight-binding Hamiltonian, and frequently works at the 
level of the hoppings only. If the system has spatially 
uniform strain, the nature of the underlying lattice is in- 
deed irrelevant for the electronic structure, and taking 
the lattice deformation explicitly into account is equiva- 
lent to taking the undeformed lattice with a rescaling of 
the momenta in the undeformed BZ. Thus, for spatially 
uniform strain it is justified to ignore the lattice deforma- 
tion, and concentrate only on the hoppings. In practice 
this entails keeping the original Si in the phase factor of 
eq. ([2]), and using the deformed Si only in the Stij. 

However, if the system is not uniformly strained, in 
general one can no longer ignore the effects of lattice de- 
formations. In a Thomas- Fermi spirit, the local (but still 
on scales safely larger than the lattice spacing) vector po- 
tential A(r) is responsible for displacing the Fermi sur- 
face in reciprocal space: E(k) — >• E(k — f A). Different 
regions experience different Fermi surface displacements, 
and one necessarily needs an absolute frame of reference 
in reciprocal space to track the position of the Fermi sur- 
face throughout the entire system. That is where the ex- 
plicit consideration of the deformed Si in the phase factor 
of eq. ([2| is important. We now look at that explicitly. 

A physically accurate parametrization of the varia- 
tion of hopping amplitude with inter-carbon distance is 
t(Si) = toexp[—f}(6i/a - 1)] with a = 1.42 A being the 



unstr ained n earest neighbor separation, to — 2.7eV, and 
(3 ~ j 14 l i7 | 3i [ rpj^ i^g^ anc [ direction of the three near- 
est neighbor vectors Si transforms under strain according 
to S[ = (I + e)'Si, where / is the 2x2 identity matrix, and 
e the two-dimensional Cartesian strain tensor, with the 
x axis along graphene's armchair direction. In reciprocal 
space, the Hamiltonian, eq. ([I]), linearized to first order 
in strain reads: 



H 
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H.c. (3) 



The new term proportional to k-e-Sj arises from expand- 
ing exp(— ik ■ S' j) to linear order in strain. Consequently, 
it contributes on equal footing with Stj, which is also of 
leading order in strain^. 

Using the undeformed BZ as global reference, and ap- 
proximating this Hamiltonian near the K points, one re- 
covers the form of the unstrained Hamiltonian but with 
the substitution (fc —> k — f A), where the vector poten- 
tials are now given by: 



A Kl = -A K > = 
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(4b) 
(4c) 
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4>o = ^, and the various Ki points are defined as in 
Fig. [ljd). The common term A p is proportional to the 
logarithmic derivative of the hopping, /?, and arises from 
the hopping perturbations, Stj, alone. It has th e ex- 
pected dependence on the strain tensor component J 14 | 29 [ 
The herein derived additional terms are the corrections 
due to lattice deformations. Since f3 ~ 3, the lattice cor- 
rections are equally important, not only for being of the 
same order in strain, but for having similar numerical 
coefficients. It is also worth emphasizing that taking ex- 
plicit consideration of the lattice deformations leads to a 
vector potential that is different for all the corners of the 
BZ. This is of course expected because under an arbitrary 
deformation the equivalence among the three K and K' 
points is lost. In fact, the expressions above quantify the 
fact that upon general deformation of the graphene lat- 
tice, the Dirac point no longer lies at Ki, nor at any of 
the symmetry points of the deformed BZ, a point empha- 
sized early on 17 . The only remaining constraint is time- 
reversal symmetry, which forces = —Ajc'.- Finally, 
it should be noted that the pseudo vector potential de- 
pends on the crystallographic orientation relative to the 
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strain and that, unlike a real vector potential, it has no 
gauge freedom as it is determined by an observable. 



III. DISCUSSION 

All existing calculations consider only the term A p and, 
as a result do not properly account for the shift in the 
Dirac points. Consider, for example, the seemingly triv- 
ial case of tensile isotropic strain. The traditional form 
of the pseudo vector potential, A p , predicts that there 
should be no shift in the Dirac points (e xx = e yy and 



in equation 4d). However, the BZ shrinks isotrop- 



ically under tensile isotropic strain, and, by symmetry, 
the Dirac points should follow the corners of the deformed 
BZ resulting in a Dirac point dependent shift toward the 
T point with respect to the unstrained reference state. 
In Fig. [2j the reciprocal space shifts of the Dirac points 
predicted by the traditional and herein corrected forms 
of the pseudo vector potential are compared. The con- 
tours are the strained band structure calculated using 
a nearest neighbor tight binding model which accounts 
for both the strain induced changes in hopping ampli- 
tudes and the lattice deformation^. For isotropic tensile 
strain, the lattice-corrected pseudo vector potential in 
eqs. Q properly predicts the displacement of each Dirac 
point due to strain. The less trivial cases of uniaxial or 
shear strain are also shown in Fig. [2] The differences 
between the red (traditional) and orange (corrected) ar- 
rows make it clear that the lattice corrections are needed 
to determine the absolute position of the Dirac point in 
reciprocal space (they can even reverse the sign of A). 

As discussed above, these corrections are immaterial 
for systems with spatially uniform strain but must be 
included for systems with non-uniform strain. For exam- 
ple, if there are regions of isotropic tension embedded in 
regions of different strain states, the relative shifts are im- 
portant, even though deformations are locally isotropic. 
In these cases lattice corrections contribute an extra, 
Dirac point specific, space dependence which is impor- 
tant when calculating transport across a strain barrier or 
when considering the spatial distribution of the pseudo- 
magnetic fields: Br. = V x 

In particular, these corrections are critical when try- 
ing to optimize the Landau quantization caused by the 
strain-induced pseudomagnetic fields. Ideally one desires 
the pseudomagnetic field to be nearly constant through- 
out the system so that the Landau levels are as narrow 
and well defined as possible. This imposes a delicate and 
non-trivial constraint on what deformations are compat- 
ible with constant pseudomagnetic fields. An early sug- 
gestion uses the traditional form of the pseudo vector 
potential to conclude that strain profiles with an over- 
all trigonal symmetry tend to generate smooth effective 
pseudomagnetic fields 7 . To show how the lattice correc- 
tion can quantitatively and qualitatively affect this con- 
clusion, we simulate a situation where graphene is cov- 
ering an equilateral triangular pit with a uniformly dis- 



tributed vertical load. This geometry was originally pro- 
posed by Guinea et al as a simple method of generating 
fairly uniform pseudomagnetic fields 7 -'. We calculate the 
strain fields using finite element analysis, and extract the 
local pseudomagnetic fields from eqs. Q. The finite ele- 
ment analysis was performed using Comsol Multiphysics, 
with a two-dimensional thin plate model including geo- 
metric non-linearity. The edges were fixed and the pres- 
sure was applied using a face load. Graphene's Young's 
Modulus of 1 TPa and thickness of 3.5 were used 
along with the Poisson ratio of graphite of 0.16£P^. To 
make the triangles more realistic we include 2 nm radius 
fillets on the three corners, which smooth down the sharp 
boundary corners. The surface was meshed with trian- 
gles with a maximum element size of 1 nm, and strain 
fields were evaluated in the mid-plane of the plate. 

The effect of the lattice corrections to the extracted 
pseudomagnetic field are shown in Fig. [3] for an equi- 
lateral triangle with 50 nm sides, and under 14 MPa of 
pressure. At this pressure the graphene sheet has less 
than 0.26 % strain. The field derived from the tradi- 
tional form of the pseudo vector potential ( |4d|_ results 
in a fairly uniform pseudomagnetic field [Fig. pFa)], in 
agreement with Guinea et aP. In contrast, Figs. |3[b,c) 
show that the new terms cause the electrons near the 
three different ^1,2,3 points to experience different pseu- 
domagnetic fields, which vary strongly across the system. 
The differences in the pseudomagnetic fields at the differ- 
ent Dirac points may be extremely useful in the context 
of valleytronics^. 



IV. CONCLUSION 

In summary, accounting for explicit lattice deforma- 
tions in the calculation of the pseudo vector potentials 
generates new leading order, and K point specific terms, 
that are needed to accurately describe the strain depen- 
dent shift of the Dirac points in reciprocal space. These 
terms are important in situations of non-uniform strain, 
such as when exploring transmission across strain bar- 
riers, or when predicting pseudomagnetic fields arising 
from particular strain profiles. In those situations the 
precise space dependence of the pseudo vector potentials, 
Ar-., and the understanding that they are different at 
each of the three inequivalent Dirac points is important. 



ACKNOWLEDGMENTS 

VMP is supported by the NRF-CRP award "Novel 2D 
materials with tailored properties: beyond graphene" (R- 
144-000-295-281). 



5 




FIG. 2. (Color online) Contours of the band structure of graphene under tensile isotropic strain, shear strain, uniaxial armchair 
strain, and zig-zag strain (rows), for e — 1% near the three K points (columns). The contours are overlaid with the Brillouin 
zone of unstrained (solid, black) and strained (dashed, brown) graphene. Vectors mark the displacement of the Dirac points 
predicted by the traditional (A p , dashed/red arrow) and the corrected (A^, solid/orange arrow) form of the pseudo vector 
potential (eqs.Q, with the green dots marking the positions of the Dirac points for strained graphene. The red vectors appears 
as a dot for isotropic strain because the traditional form of the vector potential does not predict a shift in the Dirac points. 
Each plot is square with an area of 0.12 2 . 
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FIG. 3. The spatial distribution of the pseudomagnetic fields generated when an equilateral triangle with 50 nm sides, 2 nm 
radius fillets, and the base oriented 30 degrees counterclockwise from the armchair direction is pressurized to 14 MPa. In (a) 
the calculation is done using the traditional form of the pseudo vector potential, while in (b) it is performed including the 
corrections derived in this work. Here, the three different colored curves correspond to the pseudomagnetic fields felt by the 
electrons at the three different K points with the magnetic field for the electrons at Ki isolated in (c). (Color online) 
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